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Abstract. We present an analysis of CCD photometric observations of the eclipsing novalike cataclysmic variable DW UMa 
obtained in two different luminosity states: high and intermediate. The star presents eclipses with very different depth: ~ 1.2 
\ mag in the high and ~ 3.4 mag in the intermediate state. Eclipse mapping reveals that this difference is almost entirely due 

^ . to the changes in the accretion disc radius: from ~ 0.5i?Li in the intermediate state to ~ 0.75Wli in the high state (^Li is the 

' distance from the white dwarf to the first Lagrangian point). In the intermediate state, the entire disc is eclipsed while in the high 

, state, its outer part remains visible. We also find that the central intensity of the disc is nearly the same in the two luminosity 

CO ' states and that it is the increase of the disc radius that is responsible for the final rise from the 1999/2000 low state. We find 

CO ' that the intensity profile of the disc is rather flat and suggest a possible explanation. We also discuss the effect of using a more 

>^ , realistic limb-darkening law on the disc temperatures inferred from eclipse mapping experiments. Periodogram analysis of the 

' high state data reveals "positive superhumps" with a period of 0!l455 in 2002 and 0!l461 in 2003, in accord with the results 

of Patterson et al. However, we cannot confirm the quasi-periodic oscillations reported by these authors. We obtain an updated 
^ . orbital ephemeris of DW UMa: T„,,„[HJD] = 2446229.00687(9) + O) 136606527(3)£. 

^ , Key words, accretion, accretion discs - stars: individual: DW UMa - novae, cataclysmic variables - X-ray: stars 

B 

C3 1. Introduction Biro f2000^ presented V and R eclipse maps of DW UMa in 
^ a high state. The eclipse maps show that the radial profile of the 
DW Ursae Majoris is an eclipsing nova-like (NL) cataclysmic ^^cretion disc (AD) temperature distribution is rather flat, sim- 
variable (CV) with an orbital period of ~ 3 .27 hours. The star jj^j^ ^.^^ ^^^^ ^^^^^ ^^^^^ 3^ novalikes (Rutten 
H is one of the NLs which show deep (~ 4 mag) low states ^992 Baptista et al.i996;. Biro OOQOl also estimated the 
- - . (Hessman 1990; Honeycutt et al. 1993 2003 1, and "negative" ^-^^^^^^ ^^^^^^ 270 ± 50 pc. 
and "positive superhumps" (Patterson et al. 2002 1. Shafter et 
al. (1988 1 carried out the first detailed photometric and spec- 
troscopic study of DW UMa in a high state. The photometry gpg^j^^j observations of Dhillon et al. ( fTOt in low state 
revealed ~ 1 .5 mag deep eclipses, with a peculiar V-Uke shape, revealed narrow Balmer emission lines originating from the ir- 
The spectral observations showed single-peaked H, He: and radiated face of the secondary star. Marsh & DhiUon tl9'97l 
He II emission lines with superimposed transient narrow ab- ^^^^j^ ^^^^^^ ^j^^ secondary star in their /-band spectra ob- 
components. These peculiai-ities led Thorstensen et al. ^^-^^^ ^^^^^ ^^^-^^^ ^ j^^^^ ^-^-^ jj^^ ^-^^^^^^ 450 
GM) to group DW UMa along with SW Sex, V1315 Aql and 850 pc depending on the assumed spectral type of the sec- 
PX And, and to define so-caUed SW Sex stars. ondary. Ultraviolet spectroscopy during a recent low state was 

presented by Knigge et al. (2000). Surprisingly, they found that 

Send offprint requests to: Y. Stamshe\ , ■ ^, , , , . , j ri^crv a 

. , ' , , ■ , „ , XT ■ , . ■ , during the low state the continuum shortward of 1450 A was 

based on observations obtamed at Rozhen National Astronomical ... , , , , . , , • , r^, . , , ■ , 

„, ^ Ti 1 ■ . TT u T • . ^ J . higher and bluer than in the high state. This led Knigge et al. 

Observatory, Bulgaria, at Hoher List, Germany and at Kryoneri, ° o 00 

(2000) to suggest that in the high state of DW UMa the white 
** e-mail- valKaastro bas bg (VS) dwarf (WD) is permanently hidden from our sight by the rim of 

zk(aastro. bas.bg (ZK) , henri .boffin(aoma.be (HMJB) , a flared accretion disc. Knigge et al. (2000) also modeled the 

vgenkov(aastro. bas.bg (VG) , Christina. Papadaki^oma. be ultraviolet spectrum of the WD and estimated the distance to 
(CP) , carpano@astro.uni-tuebingen.de (SC) DW UMa to be 830 + 150 pc. 
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Fig. 1. Photometric observations of DW UMa. Note the much larger range of the y-axis in the upper row panel. 



Using the same observations as Knigge et al. ('20001, 
Araujo-Betancor et al. (2003 1 could measure the eclipse phases 
of the WD. This allowed them to determine the system parame- 
ters of DW UMa with very good accuracy. Repeating the anal- 
ysis of Knigge et al. (2000 1, and using /- and /T-band photom- 
etry, Araujo-Betancor et al. (.2003f obtained two estimates of 
the distance to DW UMa: 930 + 150 pc and 590 + 100 pc, 
respectively. 

In this article we present photometric observations of DW 
UMa in two different brightness states - intermediate and high. 
We performed periodogram analysis in the high state and de- 
tected "positive superhumps". We also present eclipse maps of 



DW UMa's accretion disc in the intermediate and high states. 
In the last part of the article we discuss the results and present 
our conclusions. 



2. Observations and data reduction 

Photometric CCD observations of DW UMa were obtained 
with the 2.0-m RCC telescope at Rozhen Observatory, the 0.85- 
m Schmidt telescope at the Royal Observatory of Belgium 
(ROB), the 1-m telescope at Hoher List Observatory and the 
1.2-m telescope at Kryoneri Observatory between 2000 and 
2003. Depending on the telescope and the atmospheric con- 
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Table 1. CCD observations of DW UMa. 



Date 


Start 


duration 


mid-eclipse 




HJD-2450000 


[hours] 


[HJD] 


Roziien 2-m, Joiinson V filter 


9000 Tan 10 

Z,UWU Jail ^\J 


1 5^4 ^4 


7 1 ^ 


















1 J J^■.UlilU^■ 


9000 Tan 1 ^ 


1557.44 


4 fiO 


l S'M 480QQ 








1 S"^? 6174.^ 

1 J J / .U 1 / J 


9009 Tan 1 1 
zx)\jA, J all 11 




1 70 

1 . / V 


99Rfi 9770^; 


2002 Jan 13 


2288.32 


7.50 


2288.32641 








2288.46296 








2288.59911 


2002 Jan 14 


2289.40 


4.44 


2289.41825 








2289.55524 


2002 Feb 9 


2315.25 


9.95 


2315.37403 








2315.51110 








2315.64799 


2002 Feb 17 


2323.22 


8.86 


2323.29749 








2323.43444 








2323.57099 


ROB Schmidt telescope, unfiltered 


2000 Mar 22 


1626.38 


5.50 


1626.46748 








1626.60419 


2002 Feb 16 


2322.39 


3.96 


2322.47774 


2002 Apr 8 


2373.45 


4.24 


2373.56879 


2002 Apr 10 


2375.38 


4.43 


2375.48054 


2002 Apr 24 


2389.38 


4.63 


2389.41500 








2389.55209 


2003 Mar 16 


2715.52 


2.68 




2003 Mar 22 


2721.51 


4.23 


2721.64166 


2003 Apr 15 


2745.38 


5.00 


2745.41090 








2745.54838 


HoherList 1-m, unfiltered 


2000 Feb 10 


1585.50 


3.26 


1585.62166 


2000 Feb 1 1 


1586.60 


3.80 


1586.71482 


2001 Jan 15 


1925.37 


4.69 


1925.49886 


2001 Jan 16 


1926.58 


4.64 


1926.59172 








1926.72804 


2003 Mar 04 


2703.44 


1.52 


2703.47304 


2003 Mar 06 


2705.34 


8.31 


2705.52281 








2705.65868 


2003 Mar 07 


2706.30 


7.47 


2706.34213 


Kryoneri 1.2-m, unfiltered 


2003 Jun 02 


2793.31 


2.89 


2793.36068 



ditions, the exposure time used was between 20 and 60 s. 
A journal of the observations is given in Table[2 The CCD 
frames were proceeded in the standard way with bias removal 
and flat-field correction, followed by aperture photometry with 
the DAOPHOT procedures (Stetson 1987). The Rozhen ob- 
servations were performed with a Johnson V filter, while the 
rest are unfiltered. For Rozhen observations the stars DW 
UMa-3 (y=16.0) and DW UMa-2 (y=17.043) from Henden 



& Honeycutt ( I1995> served as a comparison and check, respec- 
tively. For the unfiltered observations as a comparison star we 
used the star GSC 3822 0070, whose standard magnitudes are 
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Fig. 2. Part of the long-term light curve of DW UMa obtained 
by RoboScope (Honeycutt et al. .2003 ) with the epochs of our 
observations marked with the vertical bars and the mean out-of- 
eclipse magnitudes shown with open squares. Our observations 
in 2000 have been obtained on the rising branch of a deep low 
state. For clarity, the RoboScope data obtained during eclipse 
have been removed. 



given by Biro (j^OOO). To be able to compare the V photometry 
with the unfiltered one, we have corrected the unfiltered mag- 
nitudes in the following way. During the observations in year 
2002 DW UMa was in a high state. Thus, we have computed 
the magnitude of GSC 3822 0070 in the instrumental system of 
the ROB Schmidt telescope so that the average out-of-eclipse 
magnitude of DW UMa in that system to be equal to the av- 
erage out-of-eclipse V magnitude computed from the Rozhen 
data. We obtained ~ 13.02 which is just between the Y - 13.29 
and R - 12.88 magnitudes estimated by Biro (2000). We fur- 
ther assumed that GSC 3822 0070 has the same magnitude in 
the system of the Hoher List 1-m telescope and corrected those 
magnitudes too. The Kryoneri photometry is relative to DW 
UMa-2. For DW UMa-2 we have obtained an unfiltered mag- 
nitude of ~ 16.36 and use it to put Kryoneri magnitudes in V. 
All the runs are shown in Fig.^ 

From Fig.[none can see that in 2001-2003 DW UMa was in 
its normal high state with out-of-eclipse magnitude V ^ 14.5, 
but in 2000 we have caught the star on the rising branch from 
a deep low state. The later is more easily seen in Fig.|2]where 
we show part of the long-term light curve of DW UMa ob- 
tained by RoboScope (Honeycutt et al. 2003 ) with the epochs 
of our observations marked with the vertical bars. In Fig.|5]the 
open squares show the nightly mean out-of-eclipse magnitude 
of DW UMa during our observations. It is seen that our obser- 



vations follow the photometry of Honeycutt et al. J2003> very 
well. This also gives us confidence in the coiTection of the un- 
filtered observations described above. 
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increase/decrease of the system brightness (22 March 2000, 16 
Jan 2001, 13 and 14 Jan 2002). The same seems to be also 



14.0 



14.5 15.0 

out-of-eclipe magnitude 



15.5 



Fig. 3. Eclipse depth as a function of the out-of-eclipse magni- 
tude of DW UMa. Filled circles - Rozhen V-band photometric 
observations; open circles - unfiltered observations. 

Figure[2 also reveals the huge difference in the eclipse 
depth in 2001-2003 and Jan 2000: ~ 1.2 and -3.4 mag, re- 
spectively. To our knowledge the eclipses we see in 2000 are 
among the deepest ever observed in a CV. The only close ex- 
ample is GS Pav, which sometimes also shows ~ 3.2-3.3 mag 
deep eclipses (Groot et al. 1998 1. Figure|31 shows the depen- 
dence of the eclipse depth on the magnitude at orbital phase 
zero. The later was estimated by fitting low-order polynomial 
to the out-of-eclipse measurements before and after the given 
eclipse. The dependence is almost linear and the eclipse depth 
decreases with system luminosity. 

3. Results 

3.1. The orbital ephemeris 

The eclipse timings listed in Table[l]were determined by fitting 
parabolas to the lower half of the eclipses. These timings were 
analyzed together with those given by Dhillon et al. (1994i, 
Biro ©00 1, Biro & Borkovits (20001 and Borkovits et al. 
(1200 1 1 ED02. .2003 1 to obtain the following orbital ephemeris: 



T^i„[HJD] = 2446229.00687(9) + O?136606527(3)£, 



(1) 



that is very close to that derived by Biro (2000). The nns around 
the fit is ~ 31 s. We note that Araujo-Betancor et al. ( 2003 ) 
used a different ephemeris, which left them with (O-C) residu- 
als of -1-73 s. We note that our ephemeris (and that of Biro lTOOOt 
nearly corrects for this residuals. 

3.2. Periodogram analysis 

After removing the measurements during the eclipses from 
the light curves, the data obtained in 2002 and 2003 were 
searched for the "positive superhumps" reported by Patterson 



et al. (I2002> . The periodogram of the "raw" light curves did 
not give any conclusive results. The most likely reason is the 
combination of very complex light curves, removing ~ 1 /3 of 
the data around ecUpse and the sparse distribution of the ob- 
servations. We would like to note one peculiarity of the DW 
UMa light curves: in some of the runs we clearly see a hnear 



present in the observations of Biro ( 120001 see his Fig. 1) too. 
This suggests that DW UMa's mean magnitude varies on time 
scale of ~ 1 - 2 days. We note that a 2.2 days (in 2002) or 2.1 
days (in 2003) periodicity could be linked to the superhump we 
have detected in the light curves. It is however still premature 
to claim that such a periodicity as been found. From our Figs.^ 
and|3] and Biro's Fig. 1, one can see that these brightness vari- 
ations are not accompanied by variations of the eclipse depth. 
This suggests a global brightening of the AD without signifi- 
cant changes of its brightness distribution or radius. 

To account for these trends we fitted straight line to all the 
runs, subtracted it from them and repeated the periodogram 
analysis. To de-convolve the computed power spectrum with 
the spectral window, the CLEAN algorithm of Roberts et al. 
(|1987) has been used. The results are presented in Fig.|4]and 
show that in 2002 data we could detect "positive superhumps" 
with a period of P^^^ - 0?1455, which is very close to that of 
Patterson et al. (2002 1. We also see two of the other three peaks 



reported by Patterson et al. (I2002> . those close to the first har- 
monics of the orbital and "positive superhumps" periods. The 
multi-sinusoidal fit to the data gives the following full ampli- 
tudes of these modulations: 0.096, 0.053 and 0.044 mag, in 
order of decreasing period. In 2003 we also detect "positive 
superhumps" with a period of P^^^ = 0^1461 but none of the 
shorter signals. The full amplitude of the modulation is ~ 0.08 
mag. 

We also searched for periodic and quasi-periodic oscilla- 
tions (QPOs) on the minute time scale. The presence of QPOs 
with characteristic periods in the range 1000^000 s in NLs is 
now a well established observational fact (for a short review 
see Patterson et al. 2002 1. Patterson et al. (2002 1 also see in 
a few SW Sex stars, including DW UMa, stable oscillations 
with periods in the same range as the QPOs. Figure|5] shows 
the mean power spectrum of DW UMa runs in double loga- 
rithmic scale. The power spectrum has a typical "red noise" 
shape with a power-law decrease of the power with frequency 
P{ f) oc y-i-59±o.02 ^jjj shows no traces of QPOs. We also do 
not see the stable oscillations with periods of 2375 s and 2974 
s reported by Patterson et al. (2002 1. The "red noise" shape 
of the power spectrum is believed to be a manifestation of the 
flickering, which is strong in DW UMa and as usual appears 
as well defined peaks with a typical duration of 5-15 min, and 
amplitudes reaching ~ 0.2-0.3 mag. The standard deviation of 
the light curves after subtraction of the best fit is ~ 0.06 mag, 
which is consistent with what is seen in other NLs (Kraicheva 
et al. ll999al [T999b Stanishev et al. l2002l . 

We also note that in some of the runs, 2 Feb 2002 and 
24 Apr 2002 for example, a hump just before the eclipse is 
clearly seen. Such humps are commonly observed in high- 
inclination dwarf novae in quiescence and are caused by the 
hot spot formed where gas stream hits the outer edge of the 
disc. The spot also makes the eclipse to appear asymmetric, 
with a well pronounced shoulder on the egress. The eclipses 
of DW UMa, including these with preceding humps are fairly 
symmetric and with no hints of hot spot. We believe that the 
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Fig. 4. Periododgram analysis of DW UMa in 2002 - left and in 2003 - right. The insets show the mean superhump shape. 



"pre-eclipse humps" seen in DW UMa are a result of super- 
humps in conjunction with strong flickering activity. 
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Fig. 5. Mean power spectrum of DW UMa in double logarith- 
mic scale. 



3.3. Eclipse mapping 

The echpse mapping method (Horne[l985j and the Maximum 
Entropy technique (Skilling & Bryan' 1984 1 was used to recover 
the surface brightness distribution of the accretion disc in DW 
UMa. The code we use (written by one of us - V.S.) allows for 
the presence of light that is never eclipsed. We assume a flat 
disc which lays in the orbital plane. Eclipse mapping assumes 
that all brightness variations result from the eclipse of the disc 
by the secondary star. Thus, any brightness variations due to 
anisotropic emission from the hot spot or superhumps cannot 
be handled by the algorithm. To account for these brightness 
variations we fitted the out-of-eclipse flux with a low-order 
polynomial function and normalized the eclipse by the fit. 



The shape of the individual eclipses is highly variable, in 
particular in the upper half of the profile. This is most likely 
due to the flickering and is not related to real changes of the 
AD structure. As mentioned before, flickering in DW UMa is 
strong and certainly can affect the upper part of the eclipse 
profile. To reduce the influence of flickering and other noise 
the eclipses were averaged in small phases bins. We note that 
even if some of the "pre-eclipse humps" discussed above are 
caused by temporary appearance of a hot spot and in these cases 
the eclipse shape variations are not caused by flickering, when 
working with an average eclipse the spot will not be detected in 
the eclipse maps. Because of the large diff'erence in the eclipse 
depth we have analysed the eclipses in the high and intermedi- 
ate state separately. In addition, the V-band and the unfiltered 
high state eclipses were also analysed separately. In the high 
state we have enough echpses to try to reduce the influence of 
the flickering even more: we averaged only the lowest 70% of 
the point in each bin. Such an approach is reasonable because 
flickering increases the observed flux. The error bars assigned 
to the mean values are the standard errors of the mean. The 
mean eclipses are shown in Fig.|SJ?. They have been scaled so 
as to have out-of-eclipse flux equal to the mean flux at orbital 
phase zero. The mean fluxes at orbital phase zero were deter- 
mined from the fits used to normalize the eclipses. 



To use the eclipse mapping method one needs to know the 
orbital inclination / and the mass ratio q - M2/M1, where Mi 
is the mass of the WD and M2 is the mass of the secondary 
star. In an eclipsing CV, / and q are related trough the duration 
of the eclipse of the WD A<p = (f>g - 0,, where 0,- and (p^ are 
the ingress and egress phases of the WD's center In NLs, (pi 
and (f>c can only be directly measured in the low state. Araujo- 
Betancor et al. ( I2003> obtained UV observations of DW UMa 
in the low state and could measure 0, and (p,,- This allowed them 
to tightly constrain the system parameters. More specificaUy 
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Fig. 6. Eclipse mapping of DW UMa. (a) mean light curves with the best fits. The dotted and the dashed lines are the fits to 
the unfiltered high state eclipses in 2002 and 2003, respectively, overplotted on the V-band high state one; (b-e) accretion disc 
eclipse maps, the WD's Roche lobe is shown with dotted line and the shadow of the secondary at orbital phases 0.0 and 0.041 
with dashed line; (f) radial intensity profiles of the accretion disc: open squares - V-band high state, filled squares - V-band 
intermediate state, and dotted and dashed lines - unfiltered high state in 2002 and 2003, respectively. 



they obtained q = 0.39 + 0.12 and / = 82° + 4°, and we used 
these values in our ecUpse mapping analysis'. 

In certain circumstances A(p can be also measured from 
high state eclipses (Wood et al.l l992> . Such measurements rely 
on the fact that in certain combinations of ; and q, at 0, and 
<^e approximately half of the disc area is eclipsed. This is the 
case for most deeply eclipsing systems, such as DW UMa, pro- 
vided that q > 0.3 and / > 80°. Thus, A(f> can be measured from 
the points of half-intensity of the eclipse. In almost all NLs 
this turns out to be the only method available to measure A(f>. 
However, the accuracy of these measurements is usually very 
low, and DW UMa itself is a good example. Published values 
of A(p measured from high state eclipses range from ~ 0.073 
(Shafter et al. 1988 ) to ~ 0.090 (Dhillon et al [T^ Biro|2000 1. 
From our high state eclipses we determine ~ 0.064 - 0.068, 
the lowest values ever measured in DW UMa. The only mea- 
surements of A(f) close to the accurate value of Araujo-Betancor 
et al. (^003 1 have been obtained when the AD contribution is 
small, ~ 0.083, by Dhillon et al. ( iil994,) in the low state. From 
the intermediate state eclipses we obtained A(p ^ 0.080, which 
is also very close to the value of Araujo-Betancor et al. ( 120031 . 

' The detection of "positive superhumps" in DW UMa indicates a 
lower value of q ~ 0.30 - 0.33. This is well in the Icr range given by 
Araujo-Betancor et al. <2003> . We tried eclipse mapping with q = 0.3, 
but this did not change the results significantly. The results presented 
below are obtained with q = 0.39 



It is clear then that measurements of A(f) from eclipses in the 
high state and system parameters determinations based on these 
have to be interpreted with extreme caution. 

The eclipses are rather symmetric, which suggests a nearly 
symmetric distribution of the AD intensity. That is why in the 
reconstructions we used a "full azimuthal streaming" default 
image (Horne J985j . The eclipses were fitted to reduced;^^ - 1 
and the fits are shown with solid lines in Fig.|^. The coiTe- 
sponding eclipse maps are shown in logarithmic gray-scale in 
Figs.l^J? - e. The outermost disc parts displayed coiTespond to 
the level where the intensity drops below 10% of the maximum, 
which we adopt as measure of the disc radius. The coiTespond- 
ing radial profiles of the disc images are shown in Fig.|6lr. The 
main results of the eclipse mapping of DW UMa can be sum- 
marized as follows: 

- the AD radius increased from ~ 0.57?l, in the intermediate 
state to ~ 0.75/?L, in the high state (/^l, is the distance from 
the WD to the first Lagrangian point Li); 

- the AD brightness distribution in both states is very flat; 

- the central intensities of the high and intermediate state 
maps are very close (Fig.|6lf). Thus, it seems that the in- 
crease of the AD radius can fully account for the long-term 
brightness variations, at least for the last one magnitude on 
the rising part of the 1999 low state. 

- in the high state a significant amount of uneclipsed light is 
detected: ~ 9% in the V-band and ~ 17% in the unfiltered 
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maps. The uneclipsed light is essentially zero in the inter- 
mediate state. 

- in the intermediate state the AD is almost totally eclipsed, 
whereas in the high state, a large part of the disc remains 
visible at mid-eclipse. This is the reason for the large dif- 
ference in the eclipse depth. 

4. Discussion 

Perhaps the most intriguing results of this study are the AD ra- 
dius variations and the nearly equal central intensities of the 
high and intermediate state maps. As high and low states are 
caused by strong variations of mass accretion rate, the later is 
difficult to explain within the standard accretion disc theory. 
However, what we formally call "intermediate state" is in fact 
a transition from the 1999/2000 low state to the high state and 
during this transition the disc might not be in equilibrium with 
the actual accretion rate. However, this is most likely not the 
case. Figurel^' shows that the recovery from the low state is 
rather slow and takes ~ 4 months, i.e. the disc has enough time 
to adjust its structure to changes of mass accretion rate. This 
suggests that at given value of the accretion rate some kind 
of saturation occurs and further increase of the accretion rate 
does not lead to increase of the emission from the inner part 
of the disc. Instead, the rise of the system luminosity is due to 
an increase of the disc radius from ~ 0.5/?Li to ~ 0.757?l,^- 
These values are in accord with the AD radii in NLs, including 
DW UMa itself, determined by Harrop-Allin & Warner ('1996) 
from high- and low-state observations. Note that the ratio of 
the area of circles with these radii is ~ 2.25. If the AD bright- 
ness distribution does not depend on radius, this will give ~ 1 
mag increase of the brightness, as observed. The reason of this 
behavior is not clear It could be that at certain value of the ac- 
cretion rate some mechanism starts removing energy from the 
inner disc, thus preventing further increase of the emission of 
this part. Thus, it appears that this issue could be tightly related 
to the mechanism responsible for the flat intensity profile of the 
disc, which we discuss next. 

The flat AD intensity profile found in SW Sex type NLs is 
still a puzzle. This problem has been extensively discussed in 
the literature, but there is no widely accepted model. Currently, 
two of the proposed models seems to be most viable. Accretion 
disc wind (or something else) can remove energy from the in- 
ner part of the disc thus causing the profile to be flatter than 
the steady-state T^ff oc r^^/'^ dependence. The second model as- 
sumes that SW Sex stars are intermediate polars and that the 
inner part of their discs is missing. In order to account for the 
other properties of the SW Sex stars it is also assumed that 
after the first impact with the AD's edge, part of the gas in 
the accretion stream continues moving above the disc and hits 
its surface again close to the WD. In the intermediate polar 
model the stream hits the WD's magnetosphere (for a review 
see HellierJOOOl. In this later model the flat profiles found in 
the eclipse mapping experiments are simply a result of the as- 

^ We stress that the fact that we see this behavior on the rising 
branch of the low state does not necessarily mean that one should ex- 
pect the opposite on the falling branch. 



sumption that the AD extends down to the WD's surface. The 
later model has recently gained additional support. Rodriguez- 
Gill et al. (2001 2001 ) reported the discovery of periodic cir- 
cular polarization in two SW Sex stars. Patterson et al. ( 2 002> 
reported possible detections of periodic signals (~ 1000 s) in 
the light curves of several S W Sex and suggested that this could 
be a manifestation of the rotation of a highly magnetic WD. 
Note, however, that in the intermediate polars the radius of the 
WD's magnetosphere or equivalently the inner radius of the 
AD Tin depends on the mass accretion rate M: r-m oc M^'^l'' . 
Thus, the higher M the lower rin. Therefore, it can be expected 
that in the high-state eclipse mapping analysis the AD bright- 
ness distribution will appear closer to the steady-state law than 
that in intermediate state. We do not see this behavior and the 
AD brightness distribution is rather flat in both state eclipse 
maps. 

Recently, Knigge et al. dMJOl found that during the 1999 
low state, the UV continuum of DW UMa shortward of 1450 A 
was higher and bluer than in the high state. The authors suggest 
that in the high state the WD and the inner hot part of the disc 
are permanently hidden from our sight by the rim of a flared 
AD. This certainly can alter the disc intensity profile of DW 
UMa and cause it to appear flat in the eclipse maps. However, 
the weakness of this model is the need of a very high AD rim; 
its height should be at least ~ 0. 15/?disc if the system inclination 
is 82°. With the present knowledge, it is difficult to understand 
how such a high rim can be formed. Moreover, it should have a 
constant height because Knigge et al. (.2000) found that in the 
low state the UV spectrum is at all orbital phases bluer and the 
ffux higher than in high state. In the classical steady-state AD 
model this would require unrealistically high mass accretion 
rate. The impact of the accretion stream with the disc can form 
a thick bulge but it could not persist along the whole disc rim. 

Here, we suggest another possibility. Bobinger et al. (I1997> 
found that the intensity profile of IP Peg's AD in decline from 
outburst is fiat. The authors explored many possible causes of 
this and concluded that the presence of an optically thick and 
geometrically ffat layer on top of the AD's surface masks its 
true intensity profile. Bobinger et al. ( I1997> suggested that this 
layer may be formed by the wind from the disc. This is even 
more relevant in NLs, which generally have hotter discs than 
dwarf novae. 

Knigge & Drew ( I1997> modeled the wind in the NL UX 
UMa and found that between the disc's photosphere and fast- 
moving part of the wind there is a transition zone of a slowly 
outffowing and relatively dense layer We suggest that such a 
layer may also be present in DW UMa, and may shade the WD 
and the inner hot part of the disc, or what is more likely, at- 
tenuate the radiation of the inner disc and the WD. Although 
this may seem somewhat speculative, it can account for (1) the 
ffat disc profile; (2) the red UV spectrum in high state and (3) 
nearly the same temperature of the inner disc in the high and 
intermediate states. However, it should be borne in mind that 
in UX UMa itself this layer does not ffatten the AD's temper- 
ature profile; UX UMa is often given as an example where the 
results of eclipse mapping are consistent with the steady-state 
Teif r^^^'^ law. The reason for the difference may be different 
parameters of the wind. 
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Fig. 7. Radial profile of the AD brightness temperature obtained by using (i) temperature dependent limb-darkening coefficient - 
symbols connected with solid lines and (ii) with constant limb-darkening coefficient of 0.6 - symbols connected with dash-dotted 
lines. The dotted curves show the effective temperature of steady-state ADs for mass transfer rates M - 10""^ - 10"^ Mq yr"' 
onto O.77M0 WD. 



Having the eclipse map calculated one can estimate the 
brightness temperature Tbr of the accretion disc assuming 
blackbody emission. This is done for the V-band eclipse maps 
by solving for Tbr the following equation: 

_ C6 J cos i Jv{A)B,i{TBRj)AdA 
^' ~ ^v{A)AdA ■ 

Here fj is the flux emitted by the y-th element of the AD, 
Tbr.] and 5j are its brightness temperature and area, respec- 
tively, B^KTbrj) is the blackbody spectral distribution, v{A) is 
the response of the V bandpass and D is the distance to the 
system. As noted by Baptista et al. (I1998t the relation between 
the brightness temperature and the effective temperature is not 
trivial and to obtain such a relation the vertical structure of 
the disc has to be known. However, the emission of the hot 
optically thick ADs in novalikes probably closely follows the 
blackbody law and Tbr should be close to the effective tem- 
perature. Optically thick ADs are however subject to an effect 
analogous to the limb-darkening in stars and the factor C takes 
this into account. This later effect greatly complicates the situ- 
ation because it is well known that 1) the real limb-darkening 
differs significantly from often used linear law and 2) C is a 
function of surface gravity, temperature and wavelength (Wade 
& Rucinski Togs', Diaz et ?l\JW^. 

To investigate the importance of limb-darkening further, we 
have calculated C by using the available synthetic stellar spec- 
tra with solar abundances and \ogg - 5 calculated by R.L. 
Kurucz^. We calculated C for V-band wavelengths, / = 82° 
and temperatures from 5000 to 50 000 K. Not surprisingly C 
turns out to depend strongly on the temperature (non-linearly) 
and changes from ~ 0.35 at T^s = 5QQQK to ~ 0.84 at 
Teff = 50 000 K. Thus, it seems that the dependence of limb- 
darkening on the temperature can affect significantly the es- 
timation of AD brightness temperature obtained from eclipse 

^ available atlhttp : // cf aku5 . harvard . edul 



mapping. This dependence would also affect the obtained ra- 
dial profile of the temperature in ADs with temperature vary- 
ing in a wide range. This would have a direct consequence on 
the inferred accretion rates and distances to the CVs studied by 
eclipse mapping. 

Although ADs structure is different from that of the normal 
stars, Diaz et al. J1996t have shown that the difference between 
limb-darkening laws obtained with Kurucz stellar spectra and 
AD model calculations is small (see Fig. 1 in Diaz et al. ll996t . 
Thus, the main source of error when calculating C is the as- 
sumption of constant surface gravity log g - 5 (log § = 5 is the 
maximal gravity in Kurucz' models). In real discs the gravity 
changes from log^ ~ 4.5 at the outer edge to log^ ~ 6.5 - 7.0 
at the inner parts. The work of Wade & Rucinski ( I1985t shows 
that in stars for given effective temperature the limb-darkening 
coefficient decreases with increasing log g, which means C in- 
creases with log^. If this holds true for ADs and \ogg > 5, 
then the temperatures one obtains with log ^ = 5 overestimate 
the real temperature of the AD, particularly in its inner part, 
where \ogg > 5; the degree of overestimation is however diffi- 
cult to estimate without model spectra with \ogg > 5. 

Since the distance to DW UMa is rather uncertain (see 
the Introduction for details) we have solved Eq.|2] for three 
distances: 300, 600 and 900 pc. Figure^ shows the results: 
symbols connected with solid lines show solutions with limb- 
darkening as determined from Kurucz' model spectra, and 
those connected with dash-dotted lines show the brightness 
temperature obtained using a linear limb-darkening law with 
a coefficient of 0.6. Also shown with dotted curves is the ex- 
pected effective temperature radial distribution of steady-state 
ADs for mass transfer rates M = 10"^° - 10"^MQyr"' onto 
Q.IIMq WD (Araujo-Betancor et al. '2003). It is clearly seen 
that the temperature of the AD in DW UMa in both high and 
intermediate state is very far from following the steady-state 
radial dependence: it is much flatter that expected. This finding 
is in accord with the results of Biro (,2000j. Also seen is that 
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the usage of more realistic limb-darkening law results in lower 
AD temperatures. The effect on the radial dependence is how- 
ever rather small, which is a direct consequence from the flat 
intensity profile. 

Although the recent studies point to a rather large distance 
to DW UMa, > 500 pc, we note that this seems inconsistent 
with the results of eclipse mapping (Biro 120001 this study). 
The main source of concern is the relatively high temperature 
(~ 20 000 K) at the outer disc regions obtained when distance 
larger than 500 pc is assumed. To achieve this, a very high mass 
accretion rate would be needed. It seems not very likely that the 
reason for this is the usage of a simplified limb-darkening law. 
The cause of the discrepancy is a puzzle. Although the depen- 
dence of limb-darkening on both the temperature and surface 
gravity is in favor of obtaining flatter radial temperature profiles 
(compared to if constant limb-darkening coefficient is used), 
this could hardly explain the very flat profile in DW UMa. The 
tests of the algorithm over simulated eclipses showed a good re- 
covery of the radial profile. Results of some of the simulations 
of flat discs are shown in Figs.|8j3-c - (a) the V-band intensity 
profile of a steady-state accretion disc with Teff oc ^-"•''s. fi^^ 
Teff oc r^'^-^T^ and (c) a combination of both. The disc elements 
are assumed to emits as a black body at given temperature. The 
emergent radiation is passed through a V-band filter, foreshort- 
ened and limb-darkened. We used a linear limb-darkening law 
with a coefficient 0.6. Gaussian noise of S/N=100 was added 
and so simulated eclipses were fitted to reduced = 1 using a 
full azimuthal streaming default image. It is seen that the recov- 
ery in the middle part of the disc is good. Figure|8]also clearly 
shows one well-known problem of the ME technique, namely 
that it tends to smear out any sharp structures in the disc. As a 
result the outer edge of the disc could not be accurately recon- 
structed and the radial profile is steeper at the outer part. On the 
other hand, in the inner disc only the very innermost pixels are 
smoothed and this is not a serious problem for the recovery of 
the radial profile. We notice that decreasing the noise level to 
S/N=300 does not help much to obtain a better reconstruction 
and the above situation does not improve significantly. 

Of more serious concern are the various systematic effects 
that may become important in the nearly edge-on systems. As 
already discussed, the mass accretion rate in NLs is probably 
not high enough for the disc rim to obscure the disc surface if 
the inclination is ~ 80°. However, if M ~ lO^^Moyr"' the 
disc rim may still be as high as ~ 0.07 - O.l/Jdisc"*- Recently, 
Smak (^002 1 suggested that in this case if the temperature of 
the rim is comparable to that of the outer disc, in the high- 
inclination NLs the rim may contribute more than 50% of the 
system light. Clearly, this would affect the eclipse profile sig- 
nificantly and when not accounted for may bias the eclipse 
mapping results. 

Smak JI994> and Bobinger et al. J 1 9971 1 have studied some 
aspects of the reliability of the reconstruction of 3D discs. Of 
particular interest is the work of Smak ( 1994i who tested the 
reliability of the reconstruction of 3D discs if a simple flat ID 
disc geometry is used in the reconstruction. Smak's work is 

This refers to the height above the disc mid-plane. The actual size 
of the rim is twice larger. 



however not based on the ME technique. As a part of the testing 
of our eclipse mapping code we have performed a similar study. 
We simulated eclipses in a binary with / = 82°, q - 0.39, 
Ml - Q.IIMq and M ~ lO^^Moyr"^ A steady-state disc, 
whose (half)thickness varies as r^^^ is assumed. The outer disc 
rim is modeled as a cylindrical surface with a temperature equal 
to that of the outermost part of the disc. The disc has a radius 
of 0.757?Li • We then simulated eclipses for three values of the 
rim height, 0.07, 0.1 and 0.13/?di,sc, which are appropriate for 
very high mass accretion rates. The V-band intensities were 
computed as in the case of the flat disc simulations^. Gaussian 
noise of S/N=100 was added and the simulated eclipses were 
fitted assuming a flat disc geometry. 

The reconstructed radial V-band intensity profiles are 
shown in Figs.|8ji-/ Generally, our results show that the radial 
profile of a 3D accretion disc may be still well reconstructed 
assuming a flat geometry, provided that the inclination is not 
very high. However, the reconstruction clearly worsens with 
increasing the rim height and in Figs. |SJ/-/ one can see the in- 
creasing contribution of the rim emission. Thus, for very high- 
inclination NLs with i > 85° this will definitely be a prob- 
lem, especially when the rim stars shielding part of the disc. 
Note that the impression for a better reconstruction of the disc's 
edge, compared to the flat disc simulations is misleading. This 
is due to the additional light from the rim, which compensate 
for the smoothing of the disc edge. For DW UMa with / ~ 82°, 
we may conclude that emission of the outer disc rim cannot be 
held responsible for the flat intensity profile, unless it is very 
bright or/and high. Note also that with increasing rim height 
the scatter of the reconstructed intensities at given radius also 
increases. The scatter in the maps of DW UMa is nearly twice 
less than that in the simulation with the rim height 0.13/?disc, 
suggesting that the contribution of the rim in DW UMa is not 
very significant. We conclude then that the flat temperature pro- 
file is not caused by the operation of the algorithm and the as- 
sumption of flat disc geometry. 

Finally, we would like to note that the SW Sex stars are 
perhaps the most complex of all types of CVs (Warner 1995}. 
The presence of overflowing stream, possibly a flared accre- 
tion disc, accretion disc wind and a magnetic white dwarf vio- 
late some of the assumptions of eclipse mapping. In addition, 
many of the SW Sex stars show superhumps, both "positive" 
and "negative". The former are now well understood to arise in 
an elliptic, progradely precessing accretion disc and will cer- 
tainly complicate the situation ever more. More generally, the 
presence of hot spots and spiral shocks should also be consid- 
ered. Detailed investigation of the effect of all these factors on 
the eclipse mapping is a very complicated task and is out of the 
scope of this article. We will only briefly discuss some of these 
topics. 

The largest contribution of the spiral shocks should be 
~ 10 - 15% (e.g. Baptista et al. EUUni When spread 



Note that in the disc geometry assumed, the angle at which given 
disc element is seen is also a function of the orbital phase. In this case 
the foreshortening and limb-darkening have to be applied simultane- 
ously and they become part of the visibility function (or projection 
matrix) applied to the disc elements. 
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Fig. 8. Reconstruction of the radial V-band intensity profile of simulated eclipses with S/N ratio of 100. Dots - the reconstruction, 
solid line - the disc profile used to simulate the eclipse, a-c - flat discs with different radial dependence of the temperature, d-f 
3D steady-state discs with emitting outer rims of different height (see text for details). In all reconstructions the eclipses were 
fitted to reduced - 1 using a flat ID disc geometry and full azimuthal streaming default images. 



over the disc, this wifl be much smaller and can hardly affect 
the radial profile. Pre-eclipse humps similar to those in quies- 
cent dwarf novae are almost never observed in SW Sex stars, 
suggesting little contribution of the hot spot, if any at all. This 
is expected though. The traces of the hot spot in the eclipses 
of dwarf novae disappear in outburst. This generally implies 
little contribution of the spot when the disc is hot. Because 
the discs in SW Sex stars are hot, the spot should little affect 
the eclipse mapping results of these objects. Probably this is 
also the case with the superhumps, whose light source in NLs 
contributes no more than 10-20% of the total system light. A 
bigger problem with superhumps would be the different disc 
geometry, for example the elliptic disc shape during the "pos- 
itive" superhumps. At least in one system, FX And (Stanishev 
et al. l2002t . the eclipse depth varies greatly with the "negative" 
superhumps phase. This effect may be due to the presence of 
precessing tilted accretion disc and would have serious impli- 
cations if the inclination is close to 90°. Thus, both superhumps 
will affect the eclipse profile. However, when averaging many 
eclipses observed at different superhump phases, the effect of 
the superhumps will be reduced and as a result average disc 
maps will be obtained. Other possible effects such as the spirals 
being thicker than the rest of the disc and any light scattered in 
the disc wind or/and the overflowing stream are diflicult to take 
into account and estimate. In conclusion, we may say that var- 
ious systematic effects may play a role in the eclipse mapping 
reconstructions of the high-inclination NLs. Thus, the results 
of such experiments should be always discussed with this in 
mind. 

5. Summary 

The results of this study can be summarized as follows: 



- Photometric observations obtained in intermediate and high 
states show very large difference in the eclipse depth: ~ 1.2 
mag in the high and ~ 3.4 mag in the intermediate state. 
Eclipse mapping reveals that this is due to variations of the 
accretion disc radius resulting in the accretion disc being 
entirely eclipsed in the intermediate state, whereas in the 
high state, a large part of the disc remains visible at mid- 
eclipse. We found that the accretion disc radius increased 
from ~ 0.5/?Li in the intermediate state to ~ 0.75/?l, in 
the high state (/?l, is the distance from the WD to the first 
Lagrangian point Li); 

- The central intensities of the high and intermediate state 
accretion disc maps are very close. Thus, the increase of 
the disc radius fully accounts for the final rise after the 
1999/2000 low state. 

- The accretion disc brightness distribution is very flat. We 
suggest that this might be due to the presence of a dense 
wind layer on top of the accretion disc surface. This layer 
shades inner hot part of the disc and the white dwarf. This 
simple model may account for: (1) the flat disc profile; (2) 
the red UV spectrum in high state, and (3) nearly the same 
temperature of the inner disc in high and intermediate state. 

- Based on synthetic stellar spectra we investigate the impor- 
tance of limb-darkening when calculating the temperature 
of optically thick accretion discs from eclipse maps. We 
find that the dependence of limb-darkening on the temper- 
ature can have a significant effect on the inferred disc tem- 
peratures, and thus should be always taken into account. 

- The reliability of the reconstruction of 3D discs if a sim- 
ple flat ID disc geometry is used in the reconstruction is 
tested. We explored the case of a steady-state disc seen at 
inclination / = 82° and possessing an emitting outer rim 
of different heights. The temperature of the rim is assumed 
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equal to that of the outermost part of the disc. It is con- 
cluded that the radial profile of such a disc may still be well 
reconstructed assuming a flat geometry, provided that the 
height of the disc rim is less than ~ O.l^disc- If the disc 
rim is larger than that, the ID reconstructed disc shows a 
larger scatter than is presently observed in DW UMa. The 
importance of other effects that may introduce systematic 
errors in the reconstruction of nearly edge-on discs is also 
discussed. 

- Periodogram analysis of the high state data reveals the pres- 
ence of "positive superhumps" with a period of P'^^ = 
O'JMSS, in accord with the results of Patterson et al. 
However, we cannot confirm the quasi-periodic oscillations 
reported by these authors. In 2003 we also detect "positive 
superhumps" with a period of P'^^ = 0?1461. 

- We obtain an updated orbital ephemeris: TroinlHJD] = 
2446229.00687(9) + O?136606527(3)£. 

In conclusion we want to emphasize the importance of do- 
ing time series photometry and spectroscopy of eclipsing SW 
Sex stars which show low states and particularly during their 
low/high state transitions. This may give us important clues 
for understanding their nature and for accretion disc physics 
in general. 
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